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I. INTRODUCTION 

Self gravitating matter shell (see 01) became an im- 
portant laboratory for testing global properties of grav- 
itational field interacting with matter. Models of a thin 
matter layer allow us to construct useful mini-superspace 
examples. Toy models of quantum gravity, based on these 
examples may give us a deeper insight into a possible fu- 
ture shape of the quantum theory of gravity (see [^). 
Especially interesting are null- like shells, carrying a self- 
gravitating light-like matter (see 0). Classical equations 
of motion of such a shell have been derived by Barrabes 
and Israel in their seminal paper . 

In the present paper we give a complete Lagrangian 
and Hamiltonian description of a physical system com- 
posed of the gravitational field interacting with a light- 
like matter shell. The paper contains two main results 
which, in our opinion, improve slightly the existing clas- 
sical theory of a null-like shell and provide an appropriate 
background for its quantized version. The first result is 
the use of fully gauge- invariant, intrinsic geometric ob- 
jects encoding physical properties of both the shell (as a 
null- like surface in spacetime — see Q) and the light-like 
matter living on the shell. We begin with a description of 
an "extrinsic curvature" of a null-like hypersurface S in 
terms of a mixed contravariant-covariant tensor density 
Q\ — an appropriate null-like analog of the ADM mo- 
mentum (cf. |7|). For a non-degenerate (time- like or 
space-like) hypersurface, the extrinsic curvature may be 
described in many equivalent ways: by tensors or tensor 
densities, both of them in the contravariant, covariant or 
mixed version. In a null-like case, the degenerate met- 
ric on S does not allow us to convert tensors into tensor 
densities and vice versa. Also, we are not allowed to rise 
covariant indices, whereas lowering the contravariant in- 
dices is not an invertible operator and leads to informa- 
tion losses. It turns out that only the mixed tensor den- 
sity Q\ has the appropriate null-like limit and enables 



us to formulate the theory of a null-like shell in a full 
analogy with the non-degenerate case. We prove Gauss- 
Codazzi equations for the extrinsic curvature described 
by this tensor density. In particular, the above notion 
of an extrinsic curvature may be applied to analyze the 
structure of non-expanding horizons (s ee P J). 

The quantity defined in Section |III| enables us to 
consider spacetimes with singular (distribution-like) cur- 
vature confined to a null-like hypersurface, and to prove 
that the Bianchi identities (understood in the sense of 
distributions) are necessary fulfilled in this case. Such 
spacetimes are a natural arena for the theory of a null- 
like matter shell. 

The second main result consists in treating the light- 
like matter in a fully dynamical (and not phenomeno- 
logical) way. All the properties of the matter are en- 
coded in a matter Lagrangian, which is an invariant 
scalar density on 5* (no invariant scalar Lagrangian ex- 
ists at all for such a matter, because conversion from 
scalar densities to scalars and vice versa is impossible!). 
The Lagrangian gives rise to a gauge-invariant energy- 
momentum tensor-density T°{,, which later — due to Ein- 
stein equations — arises as a source of gravity. Both 
Noether and Belinfante-Rosenfeld identities for the quan- 
tity r°6 are proved: they are necessary for the consis- 
tence of the theory. We stress that the contravariant 
"symmetric energy-momentum tensor" T""^ cannot be de- 
fined unambiguously, whereas the covariant tensor Tab, 
obtained by lowering the index with the help of a de- 
generate metric on S, looses partially information con- 
tained in T";,. On the contrary, the mixed contravariant- 
covariant tensor density T^b is unambiguously defined 
and contains — as in non-degenerate case — the entire dy- 
namical information about the underlying matter. 



In Section VI we use a method of variation of the total 
(gravity -f matter) Lagrangian proposed in |^ and derive 
this way Barrabes-Israel equations for gravity, together 
with the dynamical equations for the matter degrees of 
freedom. In Section VII we show how to organize the 
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gravitational and matter degrees of freedom into a con- 
strained hamiltonian system, with the ADM mass at in- 
finity playing the role of the total ("gravity -I- matter") 
Hamiltonian. Fin ally, t he structure of constraints is an- 
alyzed in Section 



VIII 



To clarify the exposition of ge- 
ometric and physical ideas some of the technical proofs 
have been shifted to the Appendix. 



II. INTRINSIC GEOMETRY OF A NULL 
HYPERSURFACE 

A null hypersurface in a Lorentzian space-time M is 
a three-dimensional submanifold 5* C Af such that the 
restriction gab of the space-time metrics g^u to S is de- 
generate. 

We shall often use adapted coordinates, where coordi- 
nate x'^ is constant on S. Space coordinates will be la- 
beled by fc, ? = 1,2,3; coordinates on S will be labeled by 
a, 6 = 0, 1, 2; finally, coordinates on St ^.= Vtf^S (where 
14 is a Cauchy surface corresponding to constant value 
of the "time-like" coordinate = t) will be labeled by 
A,B = 1,2. Space-time coordinates will be labeled by 
Greek characters a, (3, /i, v. 

The non-degeneracy of the space-time metric implies 
that the metric gab induced on S from the spacetime met- 
ric g^n has signature (0, -|-, -|-). This means that there is a 
non- vanishing null-like vector field on S, such that its 
four-dimensional embedding to M (in adapted coor- 
dinates — 0) is orthogonal to S. Hence, the covector 
Xu = Xi^'g^i, = X'^gau vanishes on vectors tangent to S 
and, therefore, the following identity holds: 



X-'gab = . 



(2.1) 



It is easy to prove (cf. [^) that integral curves of X", 
after a suitable reparameterization, are geodesic curves 
of the space-time metric 5^^. Moreover, any null hy- 
persurface S may always be embedded in a 1-parameter 
congruence of null hypersurfaces. 

We assume that topologically we have S = M} x S"^ . 
Since our considerations are purely local, we fix the ori- 
entation of the M} component and assume that null like 
vectors X describing degeneracy of the metric gab of 5* 
will be always compatible with this orientation. More- 
over, we shall always use coordinates such that the coor- 
dinate x^ increases in the direction of X , i.e., inequality 
X(x^) = X'^ > holds. In these coordinates degeneracy 
fields are of the form X — f{do — n^dA), where / > 0, 

= goA and we rise indices with the help of the two- 
,. . , . ~AB . 

dimensional matrix g , inverse to gAB- 

If by A we denote the two-dimensional volume form on 

each surface x^ = const.: 



A := del gAB , (2.2) 

then for any degeneracy field X of gab^ the following ob- 
ject 



A 



is a scalar density on S. Its definition does not depend 
upon the coordinate system (a;°) used in the above defi- 
nition. To prove this statement is sufficient to show that 
the value of vx gets multiplied by the determinant of the 
Jacobi matrix when we pass from one coordinate system 
to another. This means that vx '■— vxdx'^ A dx^ A dx^ 
is a coordinate-independent differential three-form on S. 
However, vx depends upon the choice of the field X. 

It follows immediately from the above definition that 
the following object: 

A = vx X , 

is a well defined (i.e., coordinate- independent) vector- 
density on S. Obviously, it does not depend upon any 
choice of the field X : 



A = Xido - u^Oa) 



(2.3) 



Hence, it is an intrinsic property of the internal geometry 
gab of S. The same is true for the divergence i9aA°, which 
is, therefore, an invariant, AT-independent, scalar density 
on S. Mathematically (in terms of differential forms), 
the quantity A represents the two-form: 

L := A" {da J dx° A dx^ A dx^) , 

whereas the divergence represents its exterior derivative 
(a three-from): dL := (i9aA°) dx" A dx^ A dx^. In par- 
ticular, a null surface with vanishing dL is called a non- 
expanding horizon (see [§). 

Both objects L and vx may be defined geometri- 
cally, without any use of coordinates. For this purpose 
we note that at each point x £ S, the tangent space 
TxS may be quotiented with respect to the degener- 
acy subspace spanned by X. The quotient space carries 
a non-degenerate Riemannian metric and, therefore, is 
equipped with a volume form uj (its coordinate expres- 
sion would be: u = X dx^ A dx^). The two- form L is 
equal to the pull-back of uj from the quotient space to 
TxS. The three-form vx may be defined as a product: 
Vx = a A L, where a is any one-form on 5, such that 
<X,a>= 1. 

The degenerate metric gab on S does not allow to de- 
fine via the compatibility condition Vg = 0, any natural 
connection, which could apply to generic tensor fields on 
S. Nevertheless, there is one exception: we are going to 
show that the degenerate metric defines uniquely a cer- 
tain covariant, first order differential operator which will 
be extensively used in our paper. The operator may be 
applied only to mixed (contravariant-covariant) tensor- 
density fields H"^, satisfying the following algebraic iden- 
tities: 



'H.^X" = 



H 



ab ^ 



H 



ha 



(2.4) 
(2.5) 



Vx 



X{x^) 



where Hab := gac'H'^b- definition cannot be extended 
to other tensorial fields on S. Fortunately, as will be seen. 
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extrinsic curvature of a null-like surface and the energy- 
momentum tensor of a null-like shell are described by 
tensor-densities of this type. 

The operator, which we denote by VaH^^, could be 
defined by means of the four dimensional metric con- 
nection in the ambient space-time M in the following 
way. Given H"^, take any its extension H'^" to a four- 
dimensional, symmetric tensor density, "orthogonal" to 
5, i.e. satisfying H-"-" = ("_L" denotes the component 
transversal to S). Define VqH"^ as the restriction to S 
of the four-dimensional covariant divergence V^H'^^. As 
will be seen in the sequel, ambiguities which arise when 
extending three dimensional object H"^ living on S to the 
four dimensional one, cancel finally and the result is un- 
ambiguously defined as a covector-density on S. It turns 
out, however, that this result does not depend upon the 
space-time geometry and may be defined intrinsically on 
S. This is why we first give this intrinsic definition, in 
terms of the degenerate metric. 

In case of a non-degenerate metric, the covariant diver- 
gence of a symmetric tensor H density may be calculated 
by the following formula: 



VaH"f, = daii\ - H°^rjjj 



(2.6) 



with gac,b '■= dbQac- lu case of our degenerate metric, we 
want to mimic the last formula, but here rising of indices 
of makes no sense. Nevertheless, formula (^|^) may 
be given a unique sense also in the degenerate cas e, i f 
appli ed t o a tensor density H^j, satisfying identities (2.4) 
and (2.5). Namely, we take as H"'^ any symmetric tensor 
density, which reproduces H^^ when lowering an index: 



To express directly the result in terms of the origi- 
nal tensor density H°j, we observe that it has five inde- 
pendent components and may be uniquely reconstructed 
from H"^ (2 independent components) and the symmet- 
ric two-dimensional matri x H-a b (3 i nde pendent compo- 
nents). Indeed, identities (2^) and ( |2.5|) may be rewrit- 
ten as follows: 



H 



B 
iB 



g H, 

~BC 



CB 
A 



B ' 



(2.10) 
(2.11) 

(2.12) 



The correspondence between H°f, and (H"^,Ha_b) is 
one-to-one. 

To reconstruct H"'' from H";, up to an arbitrary addi- 
tive term CX^-X^ , take the following, coordinate depen- 
dent, symmetric quantity: 

T^AB ._ -AC ~DB . ~CB „ q ~CA 

h g tlcog -n t± ^g - n tl fjg , 



(2.13) 
(2.14) 
(2.15) 



It is easy to observe that any H"'' satisfying (2.7) must 
be of the form: 



(2.16) 



The non-uniqueness in the reconstruction of H is, 
therefore, completely described by the arbitrariness in 
the choice of the value of H°°. Using these results we 
finally obtain: 



gcb 



It is easily seen, that such a tensor-density always ex- 



ists due to identities (2.4) and (2.5), but reconstruction 
of H*^^ from is not unique, because H*^^ -I- CX^'X" 
also satisfies (2.7) if H'"= does. Conversely, two such 



= 9„H«, - i (2U\ n\ - H^cl^";,) (2.17) 



symmetric tensors H'"^ satisfying (2.7) may differ only 
by CX'^X'^. This non-uniqueness does not influence the 



value of (2.6), because of the following identity impHed 
by (P: 



= {X''X'gac)i 
= X°'X'^gac,b + 2X"'gacX'^^b 



X'^X'^gacfi ■ (2. 



Hence, the following definition makes sense: 
VaH°^ := daii\ - -ti'"^gac,b ■ 



(2.9) 



The right-hand-side does not depend upon any choice 
of coordinates (i.e., transforms like a genuine covector- 
density under change of coordinates). The proof is 
straightforward and does not differ from the standard 



case of formula (2.6), when metric gab is non-degenerate. 



The operator on the right-hand side of ( 2.17 ) may thus 
be called the (three-dimensional) covariant derivative of 
H°jj on S with respect to its degenerate metric gab- We 
have just proved that it is well defined (i.e., coordinate- 
independent) for a tensor density fulfilling conditions 
( PI ) and <^%- 

Equation ( p.9[ ) suggests yet another definition of the 
covariant divergence operator. At a given point x G S 
choose any coordinate system, such that derivatives of 
the metric components gac vanish at x, i.e., gac,b(x) = 0. 
Such a coordinate system may be called inertial. The 
covariant divergence may thus be defined as a partial di- 
vergence but calculated in an inertial system: VqH"^ := 
daii-\- Ambiguities in the choice of an inertial system do 
not allow us to extend this definition to a genuine covari- 
ant derivative VcH^^. However, it may be easily checked 
that they are sufficiently mild for an unambiguous defini- 
tion of the divergence (cf. Remark at the end of Section 

0)- 
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The above two equivalent definitions of the operator V 
use only the intrinsic metric of S. We want to prove now 
that they coincide with the definition given in terms of 
the four dimensional space-time metric-connection. For 
that purpose observe, that the only non-uniqueness in 
the reconstruction of the four-dimensional tensor density 
of H'^'' is of the type CX'^X'^ . Indeed, any such recon- 
struction may be obtained from a reconstruction of H'^'^ 
by putting H'^'^ = in a coordinate system adapted to 
S (i.e., such that the coordinate remains constant on 
S). Now, calculate the four-dimensional covariant diver- 
gence Hi, := V^H^j^. Due to the geodesic character of 
integral curves of the field X, the only non-uniqueness 
which remains after this operation is of the type CX,y. 
Hence, t he restriction H;, of Hi, to S is already unique. 
Due to (2.6), it equals: 



1 



da-H.% - -n'^'gacb = VaH-^b . (2.18) 



III. EXTRINSIC GEOMETRY OF A NULL 
HYPERSURFACE. GAUSS-CODAZZI 
EQUATIONS 

To describe exterior geometry of S we begin with co- 
variant derivatives along S of the "orthogonal vector X" . 
Consider the tensor VaX^. Unlike in the non-degenerate 
case, there is no unique "normalization" of X and, there- 
fore, such an object does depend upon a choice of the 
field X. The length of X is constant (because vanishes). 
Hence, the tensor is again orthogonal to S, i.e., the com- 
ponents corresponding to fi — 3 vanish identically in 
adapted coordinates. This means that VqX'' is a purely 
three dimensional tensor living on S. For our purposes 
it is useful to use the "ADM-like" version of this object, 
defined in the following way: 



rtiX) :- -s {vx (VfcX° - S^V^X-) + S^dcA'} 



(3.1) 



where s := sgng = ±1. Due to above convention, the 
"extrinsic curvature" (5°{,(X) feels only external orienta- 
tion of S and does not feel any internal orientation of the 
field X. 

Remark: If S is a non- expanding horizon, the last 
term in the above definition vanishes. 

The last term in (3.1) is X-independent. It has been 
introduced in order to correct algebraic properties of 
the quantity vx (V^X" — (5°Vc X'^): we prove in the Ap- 
pendix ^ (see Rem ark after ( A26| )) that Q°-fj satisfies 
identities (p]j)-(2.5) and, therefore, its covariant diver- 
gence with respect to the degenerate metric gab on S is 
uniquely defined. This divergence enters into the Gauss- 
Codazzi equations which we are going to formulate now. 
Causs-Codazzi equations relate the divergence of Q with 



the transversal component of the Einstein tensor- 
density = ^Idet^l (i?^^ -(5^ii?). The transver- 
sal component of such a tensor-density is a well defined 
three-dimensional object living on 5'. In coordinate sys- 
tem adapted to S, i.e., such that the coordinate x^ is 
constant on S, we have Q\ = Q\. Due to the fact that 
Q is a tensor-density, components Q\ do not change with 
changes of the coordinate x^, provided it remains con- 
stant on S. These components describe, therefore, an 
intrinsic covector-density living on S. 

Proposition 1. The following null-like-surface version 
of the Gauss-Codazzi equation is true: 



VaQ\{X) + SVxdb 



Vx 



(3.2) 



We remind the reader that the ratio between two scalar 
densities: 9cA^ and vx, is a scalar function. Its gradi- 
ent is a CO- vector field. Finally, multiplied by the density 
Vx, it produces an intrinsic co- vector density on S. This 
proves that also the left-hand-side is a well defined, geo- 
metric object living on S. 



To prove consistency of (3_^), we must show that the 
left-hand side does not depend upon a choice of X. 
For this purpose consider another degeneracy field: fX, 
where / > is a function on S. We have: 

-sQ^ifX) = vfx i^bifX-) - 5tV,{fX-)) + 5td,K- 



= jvx [fVbX-^ + X'^dbf - SU^cX' 
= -sO''fc(X)+AV,b-(5,"A>,, , 



(3.3) 



where :— log/. It is easy to see that the tensor 

q\{ip) :- AV,, - <5,"AV,c , (3.4) 



satisfies identity (|2.4|). Moreover, g^h = —gab^'^f,, 
which proves ( |2.5|) . On the other hand, we have 



Vfxdb 



VfX 



vxdb 



Vx 



+ (9cA^)(p^, , (3.5) 



But, using formula (2.17) we immediately get: 
VaqM^) - {dcA-)^,b , 



which proves that the left hand side of (3.2) does not 
depend upon any choice of the field X. The complete 
proof of the Gauss-Codazzi equation (3.2) is given in the 
Appendix ^ . 



^ In non-degenerate case, there are four independent Gauss- 
Codazzi equations: besides Q-^, there is an additional equation 
relating Q-^^ with (external and internal) geometry of S. In de- 
generate case, vector orthogonal to S is — at the same time — 
tangent to it. Hence, Q-^j^ is a combination of quantities Q-^ and 
there are only three independent Gauss-Codazzi equations. 
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IV. BIANCHI IDENTITIES FOR SPACE-TIMES 
WITH DISTRIBUTION VALUED CURVATURE 

In this paper we consider a space-time M with distri- 
bution valued curvature tensor in the sense of Taub . 
This means that the metric tensor, although continuous, 
is not necessarily C'^-smooth across S: we assume that 
the connection coefficients F^^^ may have only step dis- 
continuities (jumps) across S. Formally, we may calcu- 
late the Riemann curvature tensor of such a spacetime, 
but derivatives of these discontinuities with respect to 
the variable produce a 5-like, singular part of R: 



sing(i?)' 



flUK 



{5l[Tl,]-6l[vU)^i^') - (4-1) 



where by 8 we denote the Dirac distribution (in order to 
distinguish it from the Kronecker symbol 5) and by [/] we 
denote the jump of a discontinuous quantity / between 
the two sides of S. Above formula is invariant under 
smooth transformations of coordinates. There is, how- 
ever, no sense to impose such a smoothness across S. In 
fact, the smoothness of spacetime is an independent con- 
dition on both sides of S. The only reasonable assump- 
tion imposed on the differentiable structure of M is that 
the metric tensor — which is smooth separately on both 
sides of S — remains continuous^ across S. Admitting co- 
ordinate transformations preserving the above condition, 
we loose a part of information contained in quantity (4.1), 
which becomes now coordinate-dependent. It turns out, 
however, that another p art, namely the Einstein tensor- 
density calculated from (^j), preserves its geometric, in- 
trinsic (i.e., coordinate- independent) meaning. In case of 
a non-degenerate geometry of S*, the following formula 
was used by many authors (see ||, ||, |l2|, p^): 



sing(g)^'^ - G^'Six^) 



(4.2) 



where the "transversal-to-S"' part of G^" vanishes iden- 
tically: 



G^" = 



(4.3) 



and the "tangent-to-5" part G"** equals to the jump of 
the ADM extrinsic curvature Q""^ of S between the two 
sides of the surface: 



(4.4) 



This quantity is a purely three-dimensional, symmetric 
tensor-density living on S. When multiplied by the one- 
dimensional density 8{x'^) in the transversal direction, it 



^ Many authors insist in relaxing this condition and assuming only 
the continuity of the three-dimensional intrinsic metric on S. We 
stress that the (apparently stronger) continuity condition for the 
four-dimensional metric does not lead to any loss of generality 
and may be treated as an additional, technical gauge imposed 
not upon the physical system but upon its mathematical param- 
eterization. We discuss thoroughly this issue in a Remark at the 
end of the present Section. 



produces th e fo ur- dimensional tensor density Q according 
to formula (^). 

Now, let us come back to the case of our degenerate 
surface S. One of th e g oals of the present paper is to 
prove, that formulae (4.2) and (4.3) remain valid also in 
this case. In particular, the latter formula means that 
the four-dimensional quantity Q^^'^ reduces in fact to an 
intrinsic, three-dimensional quantity living on S. How- 
ever, formula (4.4) cannot be true, because — as we have 
seen — there is no way to define uniquely the object Q"* 
for the degenerate metric on S. Instead, we are able to 
prove the following formula: 



G% = [Q\{X)] , 



(4.5) 



where the bracket denotes the jump of (5°j(Ar) between 
the two sides of the singular surface. Observe that this 
quantity doe s not depend upon any choice of X. Indeed, 
formula ( |3.3[ ) shows that Q changes identically on both 
sides of S when we change X and, hence, these changes 
cancel. This proves that the singular part sing(t/)°j, of 
the Einstein tensor is well defined. 

Remark: Otherwise as in the non-degenera te c ase, 
the contravariant components G"'' in formula ( |4.2D do 
not transform as a tensor-density on S. Hence, the quan- 
tity defined by these comp onen ts would be coordinate- 
dependent. According to (4^), G becomes an intrin- 
sic 3-dimensional tensor-density on S only after lowering 
an index, i.e., in the version of G°'b- This proves that 
G'^'^ may be reconstructed from G°fc up to an additive 
term CX^^X'^ only. We stress that the dynamics of the 
shell, which we discuss in the sequel, is unambiguously 
expressed in terms of the gauge- invariant, intrinsic quan- 
tity G%. 

Proofs of the above facts are given in the Appendix 
We conclude that the total Einstein tensor of our 
spacetime is a sum of the regular part'^ reg(t/) and the 
above singular part smg{Q) living on the singularity sur- 
face S. Thus 



reg(e)^ + sing(^)^, , 



(4.6) 



and the singular part is given up to an additive 
term CX^'X^5{x^). Due to (|2j), the following four- 
dimensional covariant divergence is unambiguously de- 
fined: 



1 



/j.A,c 



(4.7) 



We are going to prove that this quantity vanishes iden- 
tically. Indeed, the regular part of this divergence 
vanishes on both sides of S due to Bianchi identities: 



^ The regular part is a smooth tensor density on both sides of the 
surface S (calculated for the metric g separately) with possible 
step discontinuity across S. 
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reg ['Vf^G'^c) = 0- As a next step we observe that the sin- 
gular part is proportional to S{x^), i.e., that the Dirac 
delta contained in sing(5) will not be differentiated, when 
we ap ply the above covariant derivative to the singular 
part ( |4.2D . This is true because sing(G)'^^ — 0. Hence, 
only the covariant divergence of G along S (multiplied 
by S(x^)) remains. Another 5-like term is obtained from 
d^Q^ci when applied to the (piecewise continuous) regular 
part of Q. This way we obtain the term [Yeg{Q)-^^5{x'^). 
Finally, the total singular part of the Bianchi identities 
reads: 

sing (V^g^J = ([reg(^)^J + VaG\) 6{x') = , (4.8) 

and v ani shes identically due to the Gauss-Codazzi equa- 
tion (3^), when we calculate its jump across S. Hence, 
we have proved that the Bianchi identity Wf^G^c = holds 
universally (in the sense of distributions) for spacctimes 
with singular, light-like curvature. 

It is worthwhile to notice that the last term in defini- 



tion (3.1) of the tensor-density Q of S* is identical on its 
both sides. Hence, its jump across S vanishes identically. 
This way the singular part of the Einstein tensor density 
(4.5) reduces to: 



(4.9) 



Remark: Possibility of defining the singular Ein- 
stein tensor and its divergence via the standard for- 
mulae of Riemannian geometry (but understood in the 
sense of distribution!) simplifies considerably the math- 
ematical description of the theory. This techniques is 
based, however, on the continuity assumption for the 
four-dimensional metric. This is not a geometric or phys- 
ical condition imposed on the system, but only the co- 
ordinate (gauge) condition. Indeed, whenever the three- 
dimensional, internal metric on S is continuous, also the 
remaining four components of the total metric can be 
made continuous by a simple change of coordinates. In 
this new coordinate system we may use our techniques 
based on the theory of distributions and derive both the 
Lagrangian and the Hamiltonian version of the dynam- 
ics of the total ("gravity -I- shell") system. As will be 
seen in the sequel, the dynamics derived this way does 
not depend upon our gauge condition and is expressed 
in terms of equations which make sense also in general 
coordinates. As an example of such an equation consider 
(4.5) which — even if derived here by technique of distri- 
butions under more restrictive conditions — remains valid 
universally. We stress that even in a smooth, vacuum 
spacetime (no shell at all!) one can consider non-smooth 
coordinates, for which only the internal metric gat on a 
given surface, say {x^ = C}, is continuous, whereas the 
remaining four components gs^ may have jumps. The en- 
tire Canonical Gravity may be formulated in these coor- 
dinates. In particular, the Cauchy surfaces {x° = const.} 
would be allowed to be non-smooth here. Nobody uses 
such a formulation (even if it is fully legitime) because 
of its relative complexity: the additional gauge condition 



imposing the continuity of the whole four-dimensional 
metric makes life much easier! 



V. ENERGY-MOMENTUM TENSOR OF A 
LIGHT-LIKE MATTER. 
BELINFANTE-ROSENFELD IDENTITY 

The goal of this paper is to describe interaction be- 
tween a thin light-like matter-shell and the gravitational 
field. We derive all the properties of such a matter from 
its Lagrangian density L. It may depend upon (non- 
specified) matter fields living on a null-like surface S, 
together with their first derivatives a daZ^ and — of 
course — the (degenerate) metric tensor gat of S: 



L = L{z 



K. K 



a; 9a,h) 



(5.1) 



We assume that L is an invariant scalar density on 
S. Similarly as in the standard case of canonical field 
theory, invariance of the Lagrangian with respect to 
reparametrizations of S implies important properties of 
the theory: the Belinfante-Rosenfeld identity and the 
Noether theorem, which will be discussed in this Section. 
To get rid of some technicalities, we assume in this paper 
that the matter fields z^ are "spacetime scalars", like, 
e.g., material variables of any thermo-mechanical theory 
of continuous media (see, e.g., Q). This means that 
the Lie derivative Cyz of these fields with respect to a 
vector field Y on S coincides with the partial derivative: 

{Cyz)'' = z'^a Y'^ . 

The following Lemma characterizes Lagrangians which 
fulfill the invariance condition: 



Lemma V.l. Lagrangian density (5.1) concentrated on 
a null hypersurface S is invariant if and only if it is of 
the form: 



L = vxf{z;Cxz;g) 



(5.2) 



where X is any degeneracy field of the metric gab on S 
and /(• ] ■ ; ■) is a scalar function, homogeneous of degree 
1 with respect to its second variable. 

Proof of the Lemma and examples of invariant La- 
grangians for different light-like matter fields are given 
in Appendix |^. 

Remark: Because of the homogeneity of / with re- 
spect to Cxz, the above quantity does not depend upon 
a choice of the degeneracy field X. 

Dynamical properties of such a matter are described by 
its canonical energy-momentum tensor-density, defined 
in a standard way: 



h 



dL 



dz 



K 



(5.3) 



It is "symmetric" in the following sense: 
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Proposition 2. Canonical energy-momentum tensor- 
density T°'i, constructed from an invariant Lagrangian 
density fulfills identities (2.4) and (2.L), i.e., the follow- 
ing holds: 



T\X'' = and T.^ = T,, 



(5.4) 



Proof: For a Lagrangian density of the form ( |5.2| ) we 
have: 



VX ( ^" 



whence: 



T%gc 



'Vxfgab = Tb, 



(5.6) 



Homogeneity of / with respect to the argument (z^ d X'^) 
imphes: 



(5.7) 



□ 



In case of a non-degenerate geometry of S, one con- 
siders also the "symmetric energy-momentum tensor- 
density" r"*, defined as follows: 



ab 



dgab 



(5.8) 



In our case the degenerate metric fulfills the constraint: 
det^afc = 0. Hence, the above quantity is not uniquely 
defined. However, we may define it, but only up to 
an additive term equal to the annihilator of this con- 
straint. It is easy to see that the annihilator is of the form 
CX'^X^. Hence, ambiguity in the definition of the sym- 
metric energy-momentum tensor is precisely equal to am- 
biguity in the definition of T""^, if we want to reconstruct 
it from the well defined object . This ambiguity is 
cancelled when we lower an index. We shall prove in the 
next theorem, that for field configurations satisfying field 
equations, both the canonical and the symmetric tensors 
coincide**. This is an analog of the standard Belinfante- 
Rosenfeld identity (see [0). Moreover, Noether theorem 
(vanishing of the divergence of T) is true. We summarize 
these facts in the following: 



In our convention, energy is described by formula: H = T^o = 
px'^z^ — L > 0, analogous to H = pq — L in mechanics and 
well adapted for Hanultonian purposes. This convention differs 
from the one used in pa] , where energy is given by Too • To keep 
standard conventions for Einstein equations, we take standard 
definition of the symmetric energy- momentum tensor t°(,. This 
is why Belinfante-Rosenfeld theorem takes form r";, = —T^i,. 



Proposition 3. // L is an invariant Lagrangian and if 
the field configuration z^ satisfies Euler-Lagrange equa- 
tions derived from L: 



dL 



da 



dL 







then the following statements are true: 



(5.9) 



1. Belinfante-Rosenfeld identity: canonical energy- 
momentum tensor T°'b coincides with (minus — 
because of the convention used) symmetric energy- 
momentum tensor t""^ : 



h = gcb 

2. Noether Theorem: 

VaT\ = . 



(5.10) 



(5.11) 



Proof: Invariance of the Lagrangian with respect to 
space-time diffeomorphisms generated by a vector field 
y on 5* means that transporting the arguments (z; dz; g) 
of L along Y gives the same result as transporting di- 
rectly the value of the scalar density L on S: 



dL 



(Cyz) 



K 



dL 



dz 



K 



[Cyz) 



K 



dL 
dgac 



{CYg)ac. 



CyL . 
(5.12) 



Take for simplicity Y — (or Y°- = S^). Hence, we 
have: (Cyz)^ = z^ba = z^ ab- Applying this and rear- 
ranging terms in the above expression we obtain: 



dL 
d^ 



da 



da 



dL 

d^a 

dL 

d^a 



Z'^b 



z'^b - 5\L 



dL 

+ -T—gacfi = . (5.13) 
dgac 



Due to Euler-Lag rang ian equations (|5.9D and to the defi- 
nitions (p]3) and (5.8) of both the energy- momentum ten- 



sors, above formula reduces to the following statement: 



daT\ 



1 



ac 

T gac,b 







(5.14) 



Our proof of this formula is valid in any coordinate sys- 
tem. In particular, we may use such a system, for which 
all partial derivatives of the metric vanish at a given point 
X £ S. In this particular coordinate system we have: 

VaT\{x) = daT\{x) = . 

But WaT'^b{x) = is a coordinate- independent state- 
ment: once proved in one coordinate system, it remains 
valid in any other system. Repeating this f or all points 
X £ S separat ely, w e prove Noet her theorem ( 5.11 ) . Sub- 
tracting now (5.14) from (5.11) we obtain the following 
identity: 



rpab 
^ gab,, 



^ab ^ 
'T^ gab,, 
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which must be true in any coordinate system. Here, both 
T"'' and t"'' are defined only up to an additive term of the 
form CX"'X'', which vanishes when multipHed by gab,c- 
In the standard Riemannian or Lorentzian geometry of a 
non-degenerate metric, the derivatives gab,c niay be freely 
chosen at each point separately, which immediately im- 
plies the Belinfante-Rosenfeld identity T = — r. In our 
case, the freedom in the choice of these derivatives is re- 
stricted by the constraint. This is the only restriction. 
Hence, the Belinfante-Rosenfeld identity is true only up 
to the annihilator of t hese constraints, i.e., only in the 
form of equation ( 5.10| ). 



□ 

Remark: In non-degenerate geometry, vanishing of 
derivatives of the metric tensor at a point x uniquely 
defines a local "inertial system" at x: if two coordinate 
systems, say (cc") and (y°), fulfill this condition at x, 
then second derivatives of x° with respect to vanish 
identically at this point. Covariant derivative may thus 
be defined as a partial derivative, but calculated with 
respect to an inertial system, i.e., to any coordinate sys- 
tem of this class. In our degenerate case, vanishing of 
derivatives of the metric does not fix uniquely the iner- 
tial system. There are different coordinate systems (x") 
and (y"), for which gab,c vanishes at x, but we have: 



dx^x'' 



ix)^0 



This is why any attempt to define covariant derivative 
for an arbitrary tensor on S fails. This ambiguity is, 
however, cancelled by algebraic propertie s of our energy- 
momentum tensor, namely by identities (2^) and (2.5). 
This enables us to define unambiguously the covariant di- 
vergence of "energy-monicntuni-like" tensor-densities us- 
ing formula (2.9). 



VI. DYNAMICS OF THE TOTAL SYSTEM 
"GRAVITY + SHELL": LAGRANGIAN VERSION 

In this paper we consider dynamics of a light-like 
matter-shell discussed in the previous Section, interact- 
ing with gravitational field. We present here a method 
of derivation of the dynamical equations of the system, 
which applies also to a massive shell and follows the ideas 
of 111. 

The dynamics of the "gravity + shell" system will be 
derived from the action principle SA — 0, where 



A^AZ 



-A^Z. 



-A, 



(6.1) 



is the sum of the gravitational action and the matter 
action. Gravitational action, defined as integral of the 
Hilbert Lagrangian, splits into the regular and the sin- 
gular part, according to decomposition of the curvature: 



^«..v = T^\/R i? = -^^(reg(i?) + sing(i?)) 



Using formulae (4.2)-(4.5), we express the singular part 
of i? in terms of the singular part of the Einstein tensor: 

M sing(i?) = -sing(g) = -G^- g^,6{x'') . (6.3) 



As analyzed in Section IV, an additive, coordinate- 
dependent ambiguity CX^X'^ in the definition of is 
irrelevant, because cancelled when contracted with g^^^: 

For the matter Lagrangian i,„attor , we assume that it has 
properties discussed in the previous Section. Finally, the 
total action is the sum of three integrals: 



A = 



L' 




'I' 


■:z + 


-/ 
















Jd 




JDnS 



(6.4) 



where D is a four-dimensional region with boundary in 
spacetime M which is possibly cut by a light-like three- 
dimensional surface S (actually, because of the Dirac- 
delta factor, the second term reduce to integration over 
Dr\ S). Variation is taken with respect to the spacetime 
metric tensor g^^i, and to the matter fields living on 
S. The light-like character of the matter considered here, 
implies the light- like character of S (i.e., degeneracy of 
the induced metric: det gab — 0) as an additional con- 
straint imposed on g. 

We begin with varying the regular part L'°Jl^ of the 
gravitational action. There are many ways to calcu- 
late variation of the Hilbert Lagrangian. Here, we use 
a method proposed by one of us (see [§). It is based on 
the following, simple observation: 



= -TT-Q'"'5g^u + T7-^/\9\ g'^'SR,.. (6.5) 

iOTT iOTT 



where 



(6.6) 



It is a matter of a simple algebra, that the last term 
of (|6.5| ) is a complete divergence. Namely, the following 
formula may be checked by inspection: 



where we denote 



1 



\9\ 9'- , 



(6.7) 



(6.8) 



(6.9) 



(6.2) 



and are not independent quantities, but the Christof- 
fel symbols, i.e., combinations of the metric components 
f/pi/ and their derivatives. In the above calculations we 
use that fact that the covariant derivative Vtt of tt with 
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respect to F vanishes identically, i.e., that the following 
identity holds: 



(6.10) 



Hence, for the regular part of the curvature we obtain: 
1 „\ 1 



167r 



g\ R 



(6.11) 



We shall integrate the above equation over both parts 
and D~ of D, resulting from cutting D with the surface 
S. This way we obtain: 



(6.12) 



only by terms containing metric components and their 
derivatives along S. Jumps of these terms across S van- 
ish identically. Hence, the following is true: 

[Q^'']d{x^) = [Qn^(x^) - sing(^)^'^ . (6.18) 

Consequently, formulae ( |6.14 ), ( |6.17 ) and ( |6.3| ) imply: 

1 



^(a;^)7r 



3^„ t--^s\rl 



167r 



9^..6 sing(g)'^- 
1 



i^av + Y^sing(^)^''%, , (6.19) 



which ends the proof of (6.13). Summing up (3.12) and 
( |6.13 ) we obtain: 



ion 



(6.20) 



Now, we are going to prove that the analogous formula 
is valid also for the singular part of the gravitational La- 
grangian, i.e., that the following formula holds: 

SL:Z = -^singiOr 6g,. + sing {d^ (nr^'K'^)) ■ 

(6.13) 

To prove this formula, we calculate the singular part of 
the divergence 9^ {'^x^'^'^^^^u) ■ Because all these quan- 
tities are invariant, geometric objects {ST is a tensor!), 
we may calculate them in an arbitrary coordinate sys- 
tem. Hence, we may use our adapted coordinate system 
described in previous Sections, where coordinate is 
constant on S. This way, using (3.9), we obtain: 



smg 



where by A we denote: 



rA pK 



(6.14) 



.15) 



(Do not try to attribute any sophisticated geometric in- 
terpretation to A^j^; it is merely a combination of the 
Christoffel symbols, which arises frequently in our cal- 
culations. It has been introduced for technical reasons 
only.) The following combination of the connection coef- 
ficients will also be useful in the sequel: 



where both terms are composed of its regular and singu- 
lar part^. 

Now, we calculate the variation of the matter part 



of the action on S: 



5i„ 



-5gab 



dL„ 



-Sz 



K 



dz^a 



dgab dz^ 



-daSz 



K 



+da {pk'^Sz'') 



(6.21) 



where we used definition ( p.8[ ) and have introduced the 
momentum canonically conjugate to the matter variable 



PK 



dz^a 



(6.22) 



Finally, we obtain the following formula for the variation 
of the total ( "matter + gravity" ) Lagrangian: 



SL = - -L,eg{gr-Sg^, 
Ion 



K 



da 



9i„ 



dz^a 



- 6{x')— {G'^' - 8nT^') Sgah 

+ d4n^^''''Sr^^^) +S{x')da{pK''Sz'') . (6.23) 



It may be immediately checked that: 

1 



n'^'^SAl, = 



167r 



gtii^i 



(6.16) 



(6.17) 



In Appendix |A| we analyze in detail the structure of quan- 
tity Q. As a combination of the connection coefficients, it 
does not define any tensor density. But it differs from the 



external curvature Q{X) of S introduced in Section IH, 



In this Section we assume that both and 6z^ vanish 
in a neighborhood of the boundary dD of the spacetime 
region D (this assumption will be later relaxed, when 
deriving Hamiltonian structure of the theory). Hence, 



the above formula was proved for regular spacetimes. 
its validity was extended to spacetimes with a three- 
dimensional, non-degenerate curvature-singularity. Here, we 
have shown that it is also valid for a light-like curvature- 
singularity. 
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the last two boundary terms of the above formula van- 
ish when integrated over D. Vanishing of the variation 
6A ~ with fixed boundar y val ues implies, therefore, the 



Euler-Lagrange equations (5.9) for the matter field 



together with Einstein equations for gravitational field. 
Regular part of Einstein equations: 

must be satisfied outside of S and the singular part must 
be fulfilled on S. To avoid irrelevant ambiguities of the 
type CX"'X^, we write it in the following form, equivalent 
to the Barrabes-Israel equation: 



Sttt\ 



(6.24) 



Summing up singular and regular parts of the above 
quantities we may write the "total Einstein equations" 
in the following way: 



3^ fdL.. 



da 



5z 



K 



(6.25) 



Here, we have defined the four-dimensional energy- 
momentum tensor: T^"^ := 8{x^)t^^ with r'^" = 0. Since 
t"-'^ was defined up to an additive term CX°-X^ ^ this am- 
biguity remains and T^'' is defined up to CX^^X'^8{x^)^ 
similarly as the quantity Q'^'^ . This ambiguity is annihi- 
lated when contracted with Sg^j^u- 



VII. DYNAMICS OF THE TOTAL SYSTEM 
"GRAVITY + SHELL": HAMILTONIAN 
DESCRIPTION 

Field equations of the theory (Euler-Lagrange equa- 
tions for matter and Einstein equations — both singular 
and regular — for gravity) may thus be written in the fol- 
lowing way^: 

5L ^ d^{i:J^''^5Tl^) + 8{x'')da{pK''5z'') . (7.1) 

Indeed, field equations are equivalent to the fact that the 
volume terms (6.25) in the variation of the Lagrangian 
must vanish identically. Hence, the entire dynamics of 
the theory of the system "matter -I- gravity" is equivalent 
to the demand, that variation of the Lagrangian is equa l 
to boundary terms only. Similarly as in equation (3.14), 
we may use definition of tt-^'^'^ and express it in terms 



Formula ( [7.]| ) is analogous to formula: dL(q, q) = (pdq)' = pdq + 
pdq in mechanics, which contains both the dynamical equation: 
p = dL/dq, and the definition of the canonical momentum: p = 
dL/dq. For detailed analysis of this structure see O. 



of the contravariant density of metric t:'^". This way we 
obtain: 



flUK 



(7.2) 



Hence, field equations may be written in the following 
way: 

5L = d^{^^''5Al^)+8{x'')da{pK''5z'') . (7.3) 

As soon as we choose a (3-|-l)-decomposition of the space- 
time M, our field theory will be converted into a Hamil- 
tonian system, with the space of Cauchy data on each of 
the three-dimensional surfaces playing role of an infinite- 
dimensional phase space. Let us choose coordinate sys- 
tem adapted to this (3-|-l)-decomposition. This means 
that the time variable t ^ x^ is. constant on three- 
dimensional surfaces of this foliation. We assume that 
these surfaces are space-like. To obtain Hamiltonian for- 
mulation of our theory we shall simply integrate equation 
(7.1) (or — equivalently — (7.3)) over such a Cauchy sur- 
face Et C M and then perform Legendre transformation 
between time derivatives and corresponding momenta. 

In the present paper we consider the case of an asymp- 
totically flat spacetime and assume that also leaves of 
our (3-1-1 )-decomposition are asymptotically fiat at infin- 
ity. To keep control over 2-dimensional surface integrals 
at spatial infinity, we first consider dynamics of our "mat- 
ter -|- gravity" system in a finite world tube hi, whose 
boundary carries a non-degenerate metric of signature 
(—,+,+). At the end of our calculations, we shift the 
boundary dlA of the tube to space-infinity. We assume 
that the tube contains the surface S together with our 
light-like matter travelling over it. 

Denoting hy V := U f^ Tjt the portion o f Et which is 
contained in the tube we thus integrate ( |7.3| ) over the 
finite volume V <Z T,t and keep surface integrals on the 
boundary dV of V . They will produce the ADM mass as 
the Hamiltonian of the total "matter -|- gravity" system 
at the end of our calculations, when we pass to infinity 
with dV = Et n dU. Because our approach is geomet- 
ric and does not depend upon the choice of coordinate 
system, we may further simplify our calculations using 
coordinate x^ adapted to both S and to the boundary 
dU of the tube. We thus assume that x^ is constant on 
both these surfaces. 

Integrating ( |7.3| ) over the volume V we thus obtain: 

5 [ L = [ (tt^^SA-^^) + [ 8ix')da (Pif"<5z^) 
Jv Jv Jv 

V Jdv Jvns 

(7.4) 

where by "dot" we denote time derivative. In the above 
formula we have skipped the two-dimensional divergen- 
cies which vanish when integrated over surfaces dV and 

vns. 
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To further simplify our formalism, we denote by pk ■= 
Pk^ the time-like component of the momentum canoni- 
cally conjugate to the field variable and perform the 
Legendre transformation: 



PkSz^) ^pkSz 



K -K 



SpK + S^pKi"") ■ (7.5) 



The last term, put on the left-hand side of ( [7.4| ), meets 
the matter Lagrangian and produces the matter Hamil- 
tonian (with minus sign) , according to formula: 



PKz'' = L„ 



Pk'z'^o 



rpO _0 



T 



(7.6) 



To perform also Legendre transformation in gravitational 
degrees of freedom we follow here method proposed by 
one of us (see [^). For this purpose we first observe that, 
due to metricity of the connection F, the gravitational 
counterpart tt'^'^SA^^ of the canonical one-form pkSz^ 
reduces as follows: 



1 

16^' 



TT^'-M?.,, = -^jkiSP'' + d, ( 7r°"5 ( ^ ) ) , (7.7) 



where P^^ denotes the external curvature of E written in 
the ADM form. Similarly, the boundary term tt'^^SA-^^ ~ 
n'^'^SA^^ reduces as follows: 



^""^^l = -^a^^^Q"' + (^""^ (jSa ) ) ' (7-8) 



where Q"'' denotes the external curvature of the tube 
dU written in the ADM form. A simple proof of these 



formulae is given in Appendix D 1 



Using these results and skipping the two-dimensional 
divergencies which vanish after integration, we may 
rewrite gravitational part of (7.4) in the following way: 



{n^-5A%y+ n^-SAj^, 
V Jav 

I6n Jv 167r Jgy 



(7.9) 



The last integral may be rewritten in terms of the hy- 
perbolic angle a between surfaces S and dlA, defined as: 
a = arcsinh(g), where 



^30 



q = 



,00„33| 



(7.10) 



and the two-dimensional volume form A = \/ det gAB on 
dV, in the following way: 



1 



-»°"^)+'-"<S)-s;A*„^ (7.1) 



For the proof of this formula see Appendix D2. Hence, 
we have: 



V Jav 



(7.12) 



gabOSd 



IGtt 



av 



1- I (XSa)- 
Jav 



Now we perform the Legendre transformation both in the 
volume: 



jkl 



(gkiSP"')' - igkiSP"' - P^'Sgki + 5 [guiP 



and on the boundary: {\5a)' — {XSa — aSX) + S{Xa). 
In appendix D 3 we prove the following formula: 



(7,13) 



Then, we have 

(7.14) 

Splitting the component Q'^q of the Einstein tensor into 
regular and singular part we obtain 



(7.15) 



The regular part of Einstein tensor density reg {G^'^) van- 
ishes due to field equations. The singular part: 



(7.16) 



meets the matter hamiltonian t^q (see formula (7^)) and 
gets annihilated due to Einstein equations: 

^ / (G% - 87rr%) - . (7.17) 

Finally, we obtain the following generating formula (cf. 



ki - mSP"') + — I (XSa-aSX) 



IGtt 



av 



vns 



[p^'kSz^ - z^Sp'k) - / 9ab5Q''' 
167r Jqv 



[ {Q^^gAB~Q'"goo) 



(7.18) 
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Using results of ||] it may be easily shown that pushing 
the boundary dV to infinity and handling in a proper way 
the above three surface integrals over dV, one obtains 
in the asymptotically flat case the standard Hamiltonian 
formula for both gravitational and matter degrees of free- 
dom, with the ADM mass (given by the resulting surface 
integral at infinity) playing role of the total Hamiltonian. 
More precisely, denoting the matter momenta by 



the final formula for dV oo reads: 



(7.19) 



167r Jy 



+ / {ttkSz'' - z^'Sttk) , (7.20) 
Jv 

where H is the "total hamiltonian", equal to the 
ADM mass at spatial infinity^. 



VIII. CONSTRAINTS 



Vacuum Einstein equations outside and inside of 5* imply 
vanishing of the regular part of G'^p,- Hence, the regular 
part of the vector constraint reads: 

reg(PiV) = , 

whereas the regular part of the scalar constraint reduces 
to: 

reg (^7 R - [p^'Pki - ^P') ^) " ° ■ 

The singular part of constraints, with support on the 
intersection sphere St = VtriS", can be derived as follows. 

Singular part of three dimensional derivatives of the 
ADM momentum P^i consists of derivatives in the direc- 
tion of x^: 

sing(P,V) = sing(93P/') - Six')[Pi'] , 
so the full vector constraint has the form 

Pi''\k^[Pi^]S{x^) ■ (8.3) 



Consider Cauchy data {P^\ gki,TrK , z^) on a three- 
dimensional space-like surface Vt and denote by g'^'' the 
three-dimensional metric inverse to gki- Moreover, we 

(3) 

use the following notation: 7 :— y/det gki, R is the three- 
dimensional scalar curvature oi gki, P ■— P'^^gu and "|" 
is the three-dimensional covariant derivative with respect 
to gu ■ 

We are going to prove that these data must fulfill con- 
straints implied by Gauss-Codazzi equations for the com- 
ponents of the Einstein tensor density. Standard de- 
composition of CJ'^p into the spatial (tangent to Vt) part 
and the time- like (normal to Vt) part gives us respec- 
tively: 



-Pi\k 



.1) 



and 



2g\n'' = -7 + ( P'^'Pki - ^P"" 



l„o\ 1 

2' I 7 



Here by n we have denoted the future orthonormal vector 
to Cauchy surface Vt- 
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00 



Components of the ADM momentum P*^' are regular, 
hence singular part of the term {P^^Pki — \P^) vanishes. 
Singular part of the three-dimensional scalar curvature 
consists of derivatives in the direction of x^ of the (three- 
dimensional) connection coefficients: 

(3) 

sing(i?)=sing {d,{Tlir -VZif')) 

= 5{x^)[Tir-TZi~g''] , (8.4) 

and expression in the square brackets may be reduced to 
the following term 



.5) 



because derivatives tangent to S are continuous. But ex- 
pression in square brackets is equal to the external cur- 
vature scalar k for the two-dimensional surface St C Vt: 



^, r"p3 %kl -pm -3/1 
7 kl9 -^ml9 \ 


= -2\ 




53 ( 






= -2\ 


/~33 


dk 





7A: = ~dk 



(8.6) 



So we get 



Formula (7.2C) is analogous to formula: —dH{q,p) = pdq — qdp 
in mechanics. In a non-constrained case this formula is equiv- 
alent to the definition of the Hamiltonian vector field (p, q) via 
tiamilton equations: p = —dH/dq, and g = dH/dp. We stress, 
however, that the formula is much more general and is valid also 
for constrained systems, when the field is not unique, but given 
only "up to a gauge" . For detailed analysis of this structure see 



(3) 



sing 7 P = 2-i^/g^'i[k]5{x^) = 2[\k]6{x-') , 



and finally: 

(3) 

1 R - 



P'^'Pki - IpA - = 2[\k]5{x^) . (8.7) 
2 / 7 
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Equations (3.3) and (8.7) give a generalization (in the 
sense of distributions) of the usual vacuum constraints 
(vector and scalar respectively). 

Now, we will show how the distributional matter lo- 
cated on St determines the four surface quantities [P\] 
and [Afc], entering into the singular part of the con- 
straints. The tangent (to S) part of Q^/^i splits into the 
two-dimensional part tangent to St and the transversal 
part (along null rays). 

The tangent to St part of Einstein equations gives the 
following: 

g% = 8nS{x^)T% , (8.8) 

which, due to ( ^.l[ ) and ( ^.3|) , implies the following two 
constraints: 



(8.9) 



The remaining null tangent part of Einstein equations 
reads: 

g\X^' = 87tS{x^)t\X^' = , (8.10) 

because r^^X'' — 0. In Appendix ^ we show that this 
equation reduces to the following constraint: 




= 



We remind the reader that the singular part of t/°3 
cannot be defined in any intrinsic way. Consequent ly, we 
have only three constraints for the singular part ( 8.11 ) 
and ( |8.9| ). The fourth constraint (in a non-degenerate 
case) has been replaced here by the degene racy condi- 
tion det ^ab for the metric on S. Equations (^), (8.11) 
together with (S^) and (8/7) are the initial value con- 



straints. 
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APPENDIX A: STRUCTURE OF THE 
SINGULAR EINSTEIN TENSOR 



We rewrite the Ricci tensor: 



(Al) 



in terms o f the following combinations of Christoffel sym- 
bols (cf. ( 6.15 ) in Section 5): 



We have: 



R 



— A'^ 



(A2) 



(A3) 



Terms quadratic in A's may have only step-like dis- 
continuities. The derivatives along 5* are thus bounded 
and belong to the regular part of the Ricci tensor. The 
singular part of the Ricci tensor is obtained from the 
transversal derivatives only. In our adapted coordinate 
system, where is constant on S, we obtain: 



sing(i?^^) 



(A4) 



where by S we denote the Dirac delta-distribution and 
by square brackets we denote the jump of the value of 
the corresponding expression between the two sides of S. 
Consequently, the singular part of Einstein tensor density 
reads: 



sing(^^) := sing ( - ) = S{x')G\ 



(A5) 



(8.11) where 



:= VWl (-5^5''" - ^Si^g""^ [Alp] = [QM • (A6) 

We shall prove that the contravariant version of this 
quantity: 

sing(^)^'' = [Q^Six'') , 

is coordinate-dependent and, therefore, does not define 
any geometric object. For this purpose we are going to 
relate the coordinate-dependent quantity Q^'^ with the 
external curvature Q"fc of S. We use the form of the 
metric introduced in | [Tc| ]: 







UA 


sM + m^nA 




UA 


9AB 


rriA 




sM + m^nA 


ruA 





(A7) 



and 
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M 


s 

A/ 




7V2 * M 


.9 Ar2 


^ A/ 


" M 








* Af 






(A8) 



±1, qab is the induced Because of identity 



where M > 0, s := sgng"^ = 
two- metric on surfaces {x^ — const, — const} and 

g IS its inverse (contravariantj metric. Both g and 
gAB are used to rise and lower indices A,B ^ 1, 2 of the 
two- vectors and m^. 



Formula (A7) implies: \/\ det ^,,, ^1 — XM. Moreover, 
the object A° defined by formula (|2.3|), takes the form 



A° — XX"" where A is given by formula (2.2) and X := 
do — n'^dA- This means that we have chosen the following 
degeneracy field: X^^ — (1, ~n^,0). 

For calculational purposes it is useful to rewrite the 
two-dimensional inverse metric g in three-dimensional 

notation, putting ^ " := 0. This object satisfies the ob- 
vious identity: 



9 9cb 



= S\ - X''6\ . 



X^ab = X''X'T,ab = -X'X'^g.afi = , (A14) 

we have also labX^ — (see |l^). Now we are going to 
use the metricity condition for the connection F: 

= + 7r"''F3, = aaTT^" + ^"''A^, . (A15) 



Consequently, 



dcA" = dc (sVMs^ J - ,c - ^ab 

= -Xg^-Xb = -Xf'lab = -XI , (A16) 



Hence, the contravariant metric (A8) may be rewritten 
as follows: 

5"" = r' - - ^ i^'^X" + m^X^) , (A9) 

where := g rriB, so that m'^ := 0, and 

It may be easily checked (see, e.g., [^o), page 406) that 
covariant derivatives of the field X along S are equal to: 



VaA = -WaA - labg Oc , 



where 



and 



^A^r° 



-g{db,VaX) = ff(V„4, A) = A^F 



(AlO) 



(All) 



(A12) 



Since A is orthogonal to S, we have Aq = 0. Due to 
(A7), the only non- vanishing component of A^ is equal 
to A3 — sM . Hence, we have lab 
and, consequently. 



sMAl, 



sXlab 



(A13) 



where I = (f^lab- 



Now, we want to calculate the component = Fg^ 



^FJ^^. Because 



r{;. = aainv/|.g| = a,in(AM) , 

it is sufficient to calculate Fg^ according to the following 
formula: 



rL — 3^Tc3a — 'j^X'' i93c,a - ^3ca) 



= ^^"53c,a - A^.g°^F^,, + A^gO^Ffc,^ 
M M M 

=Wa + —m lea 

+ ^ {(A-53c),a - A_^„(g3c - mV)} 
s 1 
M M 

Finally, we obtain the following identity: 



(A17) 



Ala =Wa+Xa + J^m^k 



(A18) 



where Xa ^dahi 
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To express Q in terms of lab and Wa, we observe that: 



(A19) 
(A20) 

sg\ = . (A21) 
The missing component is much more comphcated: 



33 



iV2 



^63 



(A22) 



and depends upon 



33- 



^33 



33 



^ 3^1 



3a 



2 (.9"W+5'^W) 



(A23) 



- -93 In A + —m'^X'gab.^ ~ ^9'''933,a 
where we have used the identity 

^ff°^5afc,3 = ^3 In A . 
We are ready to prove the foUowing 

Lemma A. 2. The object Q'^b is related with Q°^b as fol- 
lows: 

sQ\ - sQ\ - hlS% + A'^xb - S\A-Xc , (A24) 



where Xc ■= 2^cln ( — 



Proof: Using (A19), (A18) and (A9) we obtain 



sQ\ =A (^ricb - ^S\lj + K'^wb - 5%K'w, 
+ A"x6 - 5\k^xc ■ 



From definition (|3.l| ) and property ( |A10| ) one can check 
that 



stj b 



= -A5%K^^ + /) + A(u;bX"+.^ l,b)+S\Xl 
= X~ricb + A''wb-d\A-w, (A26) 



so we get (A24) 



□ 



Remark: Formula ( |A26| ), together with labX'' = = 
gabX^, gives us the orthogonaUty condition Q"'bX^ = 
and symmetry of the tensor Qab '■= gacQ^b- 

Now, we would like to examine the properties of G^"^ = 
\Q^^] ■ From continuity of the metric across S we obtain 

[lab] - sM[Al,] = sM[Vl^] = X-[T,ab] = . (A27) 



On the other hand the jump of A\^ is in general non- 
vanishing. From ( [A1§| ) we have 



[Ala] = [Wa] 



Formulae (A19) - (A21) and (A27) imply: 



s[Q\]=A'^[Al,]-6\A^[Al,] 

^A'^[wb]-5\A'[w,] = s[Q\] 



Moreover, we have 



[Q\] = . 



M 



(A28) 



(A29) 



(A30) 



-sA" [[Al^] + s—X\wb] - m'[wb] 
+ (Ai\/r'' - sm"A^) [wb] . 



(A31) 



On th e other hand the jump of may be obtained from 
(PI): 



[^33] = -[53lnA]+2m''K] 
where we have used 



But 



X-[Wa] = —[X-X^gab^3]=0. 



(A32) 
(A33) 
(A34) 



Hence 



[Q°3] = sA" {-[93 In A] + m\wb]] + M\T\wb] 



(A35) 



Using these results we calculate components of \Q^ '^] — 
G^"^. From (A3C) we can easily check the property (4.3) 



(A25) G33 = [Q33] ^ ^33 [q3^] ^ ^3fc[Q3^] ^ q 



where we used the property [Q^fc] — [Q°'b] which is crucial 
to admit that the object G°b is a well defined geometric 
object on S. On the contrary, the object G°^ is not a 
geometric object because depends on a choice of coordi- 
nates. This can be seen when we calculate the component 
G°0: 

G0"=[Q™]=.9"'[Q°3]+5°''[QM 
= ^(_[a3lnA]+m''M) 



[53A] . 



(A36) 
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It may be easily checked (see that the above quantity 
transforms in a homogeneous way with respect to coordi- 
nate transformation on S. This proves that the compo- 
nents G"'' do not define any tensor density on S. An in- 
dependent argument for this statement may be produced 
as fohows. Begin with a coordinate system in which we 
have X = do (i.e., = 0) and perform the following 
coordinate transformation: 

x° ^x° + hAX^ , = a;-^ , =x\ (A37) 

where are constant. According to (Rq) we have: 



g{dx",dx'^) = g{dx^, dx^) + bAg{dx^, dx^) 



M 



— (l-bAn^) - — , 



(A38) 



whence we get M = M. Moreover, the new tetrad 
{X, dg,di) may be calculated as follows: 



X ^X , 

dx^ ^ dx^ 

— -9o H ■■ 

dx^ dx^ 



(A39) 



= + ^Ba = S^dA - bgX , (A40) 



= 83 



This implies A = A , and, consequently, 



(A41) 



On the other hand, we have da;" — dx^ + 5^da;'^ and 



det 



\as'' J 



— 1. Hence, 



qOO _ qOO ^ G(d5°, d£°) - G(da;", da;°) 
= 26aG0^ -I- G^^bAbB , 
which does not need to vanish in a generic case. 

APPENDIX B: GAUSS-CODAZZI EQUATIONS 

We begin with the Lie derivative of a connection T with 
respect to a vector field W (see jlTt): 



(Bl) 



For the coordinate field W = da (i.e., = (5^), Lie 
derivative reduces to the partial derivative: Cw^uf = 



da . Hence, taking appropriate traces of ( |Bl| ) and 



denoting tt'^'^ :— y'\g\g'^'^ we obtain: 

^'-a,A^, = (<5>^''-<5^"''Rr^, 



We apply this formula for a = 3. This way we have: 

= db[ ( V'W^^ - VF'') } + 27^\ . (B2) 
where Tl^a := ^/\g\R^a■ But 



^ ^i** ^ afi 



Hence: 

1 



= 2g''ML + g-r;^^ + 5''^, , 

and, consequently, 

^/\g\{V^W^ - V^M/*) = 27r''-^Ai„ -I- 7r3% 
Inserting this to (p32|) we obtain: 



(B3) 



But 



2~ .a^Mi' 



(B4) 



{9''-daV\9\ + V\9\ 9^"9''^9c.p,a) 



(B5) 



where we used definition ( 6.16| )), namely 



(B6) 



Hence, we obtain the following identity: 



g^'a + db \ Q\ + ls\n'^A - ^Q^^gap^a = . (B7) 



To calculate the last term of (B7) we use the following 
Lemma B.3. The following equality holds 

1. 



sQ''^go.0,a ^X{g''g'%d - ^l9n9bc.a 



''9'\d - ^Ig'n 

^ ^cd I \ c ^bd \ d ^cb\ a3 



+ (A'g-^ + A-g"" ~ A^g-'')Ai,gtc,a 



+ 



2sQ% (da 



In M+—mBn%] . (B8) 
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Proof: From ( |A20D and ( [A2lD we obtain 
Q^^ 

and 



5 A )Wc 



where we have used formula 

Formula ( Bll ) is equivalent to ( 3^ ) if we use ( |aT^ ), and 
keep in mind the "gauge" condition X{x'^) — 1, used 
thoroughly in this proof. 



so 



gal3,a = '2Q^39^''g3b,a + Q'"'gbc,a ■ 



be. 



Moreover, from (A7) - (A8) we have 



g-^'g^bM = da In A/ + j-^men'^a 



and 



y — (,0 c5 +5 5 gscjQ d 



Af2 



33 • 



Using ( A19|) and taking into accoimt that X'^X^gab,c — 
we get 



sQ'^gbca ^Kg^'g'Xd - ^igngbc.a 



(B9) 



and finally 

sQ°''^gaf3,a =Kg^'' g'"^^ed - -jig^'')gbc,a 

+ (AV'^ + Ay* - AV')^L56c,a 



2sQS {da 



lnAI+—mBn%j . (BIO) 

□ 



Now, the proof of ( p. 2D is roughly a straightforward 
calculation starting from equation (B7) and consequent 
reexpressing all ingredients in terms of the connection 
objects lab, Wa and the metric objects M, m"*, N, X"-, gab 
describing the 4-dimensional metric (/^i,. It turns out that 
t he t erms c onta ining M, N and d rop o ut. Ins erting 
( [A2| ) and (g) into @ and using ( |A18| ), ( |A20| ), and 
(A3), we obtain: 



S^Q°''^gal3,a 



=db {-sQ^ + - K'^Xb + <5%A^Xc} 

-iA/(a,lnAf +^mBn^J 



+ ^gbc,a (A^g^'* + A^5'"^ - A'^s^^^ 

/ * R, 

X y^d + Xd + J^ril iBd 

+ 2^gbc,a (^g'"'g'"'led - ■^ig'"' 

= - sdbQ'a + ^sQ'-gbca + Xdal , 



(Bll) 



APPENDIX C: PROOF OF LEMMA ( |vTl| ) AND 
EXAMPLES OF INVARIANT LAGRANGIANS 



Since the matter Lagrangian (5.1) is an invariant scalar 
density, its value may be calculated in any coordinate sys- 
tem. For purposes of the proof let us restrict ourselves to 
local coordinate systems (x") on S which are compatible 
with the degeneracy of the metric, i.e., such that X := do 
is null-like. 

Suppose that (a;°) and (?/°) are two such local systems 
in a neighborhood of a point x G S. Suppose, moreover, 
that both vectors do coincide. It is easy to see that these 
conditions imply the following form of the transformation 
between the two systems: 



y 



y 



(CI) 
(C2) 



Three-dimensional Jacobian of such a transformation is 
equal to the two-dimensional one: det{dy^ / dx^) . Ob- 
serve that the two-dimensional part gAB of the metric gab 
transforms according to the same two-dimensional matrix 
and, whence, its determinant A gets multiplied by the 
same two-dimensional Jacobian when transformed from 
(x") to (2/°). So does also the volume vx- This means 
that the function 



/ 



L 

Vx 



(C3) 



does not change value during such a transformation. A 
priori, we could have: 



/ = /(z^;z^o;^^A;5ab) , 



(C4) 



but we are going to prove that, in fact, it cannot depend 
upon derivatives a- For this purpose consider new 
coordinates: 



y 

y° 



eix^ — £23;^ 



(C5) 
(C6) 



This implies that 



d 



d 



dy^ dx^ 



dx^ 



Passing from (x") to (y"), the value of a will be thus 
replaced by z^ a + f- AZ^ 0, whereas the remaining vari- 
ables of the function ( |C4D (and also its value) will remain 
unchanged. This implies the following identity: 



/(z^ ; ; a; gab) = /(z^ ; 0; a + eA^"" 0; gab) 



(C7) 
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which must be vahd for any configuration of the field 
. Such a function cannot depend upon a\ But in 
our coordinate system we have z^o = z^ aX'^ = Cxz^ ■ 
Thus, we have proved that 



f ^ fiz'^-Cxz'^^gab) ■ 



(C8) 



Relaxing condition ( |C2| ) and admitting arbitrary time co- 
ordinates we easily see that the dependence of (C8) 
upon its second variable must annihilate the (homoge- 
neous of degree minus one) dependence of the density vx 
upon the field X in formula ( |5.2D . This proves that / 
must be homogeneous of degree one in Cxz^ . 

As an example of an invariant Lagrangian consider a 
theory of a light-like "elastic media" described by mate- 
rial variables z^ , A — \,2^ considered as coordinates in a 
two-dimensional material space Z, equipped with a Rie- 
mannian "material metric" ^ab ■ Moreover, take a scalar 
field ^. Then for numbers a and /3 > 0, satisfying iden- 
tity 2a + (3 = 1, and for any function ip of one variable, 
the following Lagrangian density: 



and g'^' is the 3-dimensional inverse with respect to the 
induced metric gki on V. 

Let us exchange now the role of and . Identities 



(Dl) and (D2) become constraints on the boundary of 
the world-tube dU: 



1 



Al, = {dan'- + Al,n^') , 



A' 



27r33 



They imply: 



1 



167r 

where we have denoted 



gabSQ""' + da TT^'S 



_3a 



^33 



2"^^= V|detff,d|(Lr'-i"') , Lab 



(D5) 
(D6) 



, (D7) 



^33 



lb 1 



(D8) 



(C9) 



fulfills properties listed in Lemma (V.l) and, therefore, is 
invariant. If tp is constant, a possible physical interpreta- 
tion of the variable ^ as a "thermodynamical potential" , 
may be found in [ [l4| . 



APPENDIX D: REDUCTION OF THE 
GENERATING FORMULA 

1. Proof of formulae ( |7.7| ) and ( [r.8| ) 

We reduce the generating formula with respect to con- 
straints implied by identities VfeTr*"" = and Vfe7r°" — 0. 
In fact, expressing the left-hand sides in terms of tt'^'^ and 
■A^Ilv immediately get the following constraints: 



A' 

JT-l 



1 



A"o,^^^{d,n"^ + 2Al,n"^) 



(Dl) 
(D2) 



It is easy to see that they imply the following formula: 



1 



-^<?,i<5P'='+afe(^""<5(-^) ) , (D3) 



where we have denoted 



^Ofe 



(D4) 



kl 



and g""^ is the 3-dimensional inverse with respect to the 
induced metric gab on the world-tube. 



2. Proof of formula |7.1l| 



Write the right hand side as follows: 



^30 



1^00 y \^^33y VI I ^1^00^33 1 



(D9) 



and 



2ypo^^ J_^^-^^ 1 V^tg^ 
IGtt 87r ^\ + g2 



(DIO) 



This automatically implies 



33. /^^'°\ _ ^ '^g _ ^ 

^00 \^3z) g^^YT^ 87r' 



(DU) 



Proof of formula 



7.13| 



To prove ( 7.13 ), consider first the following identity: 

O/c ' 



guP'^' ^ I D + I dk[ TT^'do 



IV JV 

where we denote: 



T-OO 



(D12) 



= Col = i-xl , 
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with X = gin-, i.e., = 5q and Cx being the Lie 
derivative with respect to the field X: 

(due to Bianchi identities the right hand side is automat- 
icahy symmetric with respect to lower indices). Hence 

D : = {S'.TT^'^ S^yni'^.W.X' X'^R' 

{V^(V^X" - V°X^) + 2R"^X''} = 



16n 

(D13) 



The covariant derivative has been replaced in the 
last equation by the partial derivative 9^, because they 
both coincide when acting on antisymmetric, covariant 
bivector densities. We use also identity 



which finally implies: 



(D14) 



(D15) 



D is regular, because singular expressions contained in 
its definition cancel out, as implied by equation ( D12 ). 
Hence, we treat D as a regular expression, and there is 
no need to integrate it in distributional sense. Hence we 
have: 



gkiP"' 



8vr jy 



\g\R% 



1 

16^ 



dV 



g\ [ g'^r^.y^rf 



TT do 



^03 



^00 



From definition of a we also have: 



Ad = Stt ( 7r°°ao 



TT 

Too 



^30 



^ TT"" J \ TT 

Using the above formula we may write 



33 



(D16) 



(D17) 



i-/ Ad=i- 

JdV STTjy 



g\R\ 

30 



(D18) 

The left-hand side of the above equation is regular, but 
on the right-hand side singular terms like Jy \/\g\R^o 
and w/ay vIsKs'^rO^-r/O'^rSJ arise." The latter 
quantity, although it is a boundary term, origins from 

the volume term ^ a. (/^ (g-/^rg^ - g^^^T^^) 



Stokes theorem. From derivatives in x^-direction there 
come singular terms, which cancel out the singular part 
of i^^o, giving regular expression as a final result. 

We may rewrite expres sion s in (D18) in terms of the 
quantity Q"'' (defined by (|d|)) 



1 

16^ 



f {Q^^'gAB ~ Q°"5oo) (D19) 
167r Jqv 



what completes the proof of formula ( 7.13 ) 



APPENDIX E: PROOF OF THE CONSTRAINT 



(8.11) 



Using the decomposition (A7), (A8) of the metric, one 
can express vector n orthonormal to Vt as follows: 



If , N'^ A N"^ 



Choosing X ^ do — n^dA, we have: 



1 N 

-X ^s — {ds-m^dA) + n 



(El) 



Consequently, we can rewrite the left-hand side of ( 8.10| ) 
as follows: 

Expressing CJ'^p in t erm s of t he c anonical AD M m omen- 
tum Pki (equations ( ^.l|) and ( |8.2| )), equation ( S.IO ) takes 
the form: 



(E3) 



Equations (Q and give us the following result 

N 



APa""]) + [Afc] - 



(E4) 



Due to (A7), one can express the three-dimensional in- 
verse metric g*^' as follows: 



g'^' = (- 

^ ^ M 





A \ 

f m niA j g 


—m 




1 



(E5) 
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The above form of 5^^' can be used to rewrite the canon- 
ical momentum part of (E4): 



and finally we obtain the constraint (8.11 



N 



331 



P 



.33 



/^33 



(E6) 
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